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Abstract 

Tlie fluorescence technique has been successfully used to detect ultrahigh energy cosmic rays 
by indirect measurements. The underlying idea is that the number of charged particles in the 
atmospheric shower, i.e, its longitudinal profile, can be extracted from the amount of emitted 
nitrogen fluorescence light. However the influence of shower fluctuations and the very possible 
presence of different nuclear species in the primary cosmic ray spectrum makes the estimate of the 
shower energy from the fluorescence data analysis a difficult task. We investigate the potential of 
shower size at maximum depth as estimator of shower energy. The detection of the fluorescence 
light is simulated in detail and the reconstruction biases are carefully analyzed. We extend our 
calculations to both HiRes and EUSO experiments. This approach has shown some advantages to 
the reconstruction of the energy when compared to the standard analysis procedure. 

PACS numbers: 96.40.Pq 



* Electronic address: | vitor @ astro . lag . usp .br| 
Electronic address: "gnitanco@gmail.coni| 
t Electronic address: jortiz@astro.iag.usp.'br| 
^Electronic address: lfederico.sanchez@mi.infnltl 



I. INTRODUCTION 



Cosmic rays (CR) are the highest energy particles present in nature with energies exceed- 
ing 10^°eV . Their origin, nature and possible acceleration mechanisms are still a mystery, 
despite the efforts of many cosmic ray experiments in the last decades. Such challenge is in 
part due to the very low flux of high energy and the consequent few observed events to be 
analyzed. 

Measuring extensive air showers (EAS) is currently the only way to study the cosmic 
ray spectrum and chemical composition at energies around and above PeV. At energies 
£'>10^'' eV the shower development can be directly observed by measuring the fluorescence 
light from atmospheric nitrogen that is excited by the ionization of the secondary charged 
shower particles (essentially electrons and positrons). Experiments applying this technique 
can determine the depth of maximum air shower development (Xmax) and the corresponding 
number of chargedparticles (A^max)- Presently, the HiRes pj, Pierre Auger Observatory 
EUSO 0], OWL y and Telescope Array p]) experiments are using or planning to use 
fluorescence detectors to investigate the ultra high energy cosmic rays. 

The total amount of emitted fluorescence light in a shower is a very good approximation 
to the total number of charged particles N{X), where X is the atmospheric depth. In this 
sense the number of particles at shower maximum can serve as an estimator of the shower 
energy. The total energy that goes into electromagnetic charged particles is obtained by 
integration of the shower longitudinal profile 



where a is the average ionization loss rate pi], and the integral on the right-hand side 
represents the total track length of all charged particles in the shower projected onto the 
shower axis. 

As an alternative proposal [3| the electromagnetic energy can also be calculated by using 
the fluorescence light intensity and the fluorescence efficiency, without the need of recon- 
structing the number of particles as a function of the atmospheric depth. Such approach is 
taken as a very precise measurement of the primary shower energy because it is supposed to 
be weakly dependent on the simulation models and on the primary particle type. However, 
when the shower development details are taken into account the calorimetric measurement 
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can lead to high systematic uncertainties. Of no less concern is the important fact that the 
fluorescence efficiency as a function of air pressure, density and humidity is only known to 
a certain extent. On approach given by Eq. ^ the average ionization loss rate is used in the 
air shower reconstruction and hence the energy spectrum of the electron shower particles 
must be known via Monte Carlo simulation. 

Although the electrons and positrons constitute the majority among the charged particles 
in a shower and contribute most to the production of fluorescence light, an also important 
fraction of the shower energy is carried by particles which cannot be measured by the fluo- 
rescence technique, i.e., particles that are invisible to fluorescence telescopes. This so called 
"missing energy" has been calculated by Monte Carlo air shower simulation and contributes 
to the uncertainties involved in the method, being sensitive to primary composition. 

Theoretical works have shown the existence of a clear relation between the primary energy 
and the maximum number of particles in the shower. Recently, Alvarez-Mufiiz et al. [8| 
have studied the A^max shower quantity as an estimator for the primary shower energy, 
confirming the efficiency of this technique. Such approach was analyzed for different primary 
particles and energies using a fast one-dimensional simulation program. However, telescopes 
particularities and reconstruction procedures must be considered due to the introduction of 
biases and fluctuations in the calculation of A^max- 

The scope of this work is to explore the possibility of estimating the primary shower energy 
based on A'max, taking clS CclSG studies the HiRes [jj and EUSO j2| fluorescence experiments, 
i.e., ground and space based experiments, respectively. The telescopes particularities and 
the reconstruction procedures are included in our analysis configuring a very realistic context 
for the application of the technique. 



II. TVmax FLUCTUATIONS 

Fluctuations are intrinsic to any extensive air shower and are a cause of uncertainty in the 
energy reconstruction. In addition to that, the energy of the shower is calculated without 
the knowledge of the mass of the primary particle which initiated the shower. The nature of 
the primary particle affects the longitudinal development of showers and can exert influence 
on the reconstruction of the air shower energy, especially when primary photons are taken 
into account. 
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On the other hand, reconstruction procedures must relay on simulation programs in order 
to relate the measurable parameters to energy. In addition to the choice of the particular 
simulation program to be used, there is a general agreement that most shower reconstruction 
uncertainties originate on uncertainties on the hadronic interaction models. 

The usual procedure used by the HiRes and Auger collaborations to reconstruct the 
shower energy with fluorescence telescopes correlates the integral of the energy deposited in 
the atmosphere to the total shower energy. As mentioned before, a certain amount of the 
shower energy is carried away by particles which are invisible to the fluorescence technique, 
i.e., muons, neutrinos and high energy hadrons which are not converted to fluorescence 
photons. Such "missing energy" has been estimated by Monte Carlo simulations and shown 
to be dependent on the hadronic interaction model, primary composition, shower energy 
and arrival direction. 

In reference jsl. Song et al. have estimated that the fluctuations due to the type of primary 
particle are about ~ 5% for showers initiated by proton and iron nuclei and around ~ 20% 
if primary photons are taken into account. Meanwhile Barbosa et al. ^ and Alvarez-Muniz 
et al. have calculated similar values for hadronic primaries (~5%) which decrease with 
energy. According to 0, the hadronic interaction model is the main source of systematic 
fluctuation in the estimation of the missing energy at the highest energies. Due to the 
differences in the ■„u.«p.,ci«es of secondary part.cles sinrnlated by QGSJETOl Q and 
SIBYLL2.1 PJI hadronic interaction models, the unseen energy calculated with QGSJET is 
about 50% higher than the value predicted by SIBYLL at 10^° eV, which can be translated 
into an uncertainty of ~ 5% in the shower energy reconstruction. 

The dependence of the missing energy on the energy of the primary is a consensus among 
the various studies made so far. In addition, Barbosa and co-workers jol have studied the 
dependence of the missing energy on the arrival direction of the shower and limited it to be 
at most 0.7%. 

Furthermore, it has been suggested in reference Q that the discrepancies among such 
studies are on the order of 1-3% and could be explained by different hadronic interaction 
models for lower energy particles. It might be relevant to mention that both groups used 
different simulation programs: Song et al. and Barbosa et al. used different versions of the 
well tested CORSIK A jl2l| . while Alvarez-Muniz et al. used a hybrid fast one-dimensional 
simulation program |l3l |. 
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Considering all the uncertainties in the determination of the missing shower energy listed 
above, one can estimate that the total systematic fluctuation related to the primary energy 
reconstruction is about 5% to 10%. 

In order to investigate the fluctuations in the maximum number of particles due to the 
natural fluctuations of the shower, the different primary particles and the artiflcial fluctua- 
tions introduced by different hadronic interaction models, we have used the recently released 
CONEX shower simulation program. We have simulated 5,000 showers initiated by 
proton, iron and gamma primaries at each considered energy and for the QGSJETOl and 
S1BYLL2.1 interaction hadronic models. Shower have been simulated with the minimum 
cuts available, 1 GeV for hadrons and 1 MeV for electromagnetic particles. 

Fig. n shows the A^max distribution for 5,000 showers initiated by protons, iron nuclei 
and gamma at 10^^'^ eV simulated with the QGSJET hadronic model. The distribution 
illustrates the previous discussion regarding the fluctuations due to the primary particle 
type. If only proton and iron nuclei are considered, the iVmax distribution shows a very 
narrow distribution with median value of 1.90 x 10^° and a dispersion of 3.6% at the 68% of 
confidence level. If gamma shower are considered, the median of the combined distribution 
(gamma+proton+iron nuclei) increases slightly to 1.92 x 10^° while the dispersion reaches 
7% at 68% of confidence level. 

The same calculation has been done for proton and iron nuclei at 10^°'^ eV and the 
fiuctuations presented the same level of 3.6%. Gamma air showers have not been studied 
at the energy of 10^°'^ eV due to the fact that the CONEX program does not include 
the pre-shower algorithm effect 3, 1^ that takes into account photon interactions with 
the geomagnetic field affecting the longitudinal development of showers with energy above 
lO^^-^ eV. 

In order to investigate the uncertainty due to the hadronic interaction model we have 
simulated 5,000 proton showers with the CONEX program using SIBYLL and QGSJET 
hadronic interaction models. Fig. |21 illustrates the A^max distribution for proton showers 
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simulated with QGSJET and SIBYLL at 10 eV and the corresponding fiuctuation is 4. 
Fig. 12 also shows the same distribution for iron nuclei at 10^°'^ eV and the corresponding 
fiuctuation is 3.3%. 

Comparing the numbers given above and Figs. ^ and |21 we are able to conclude that the 
systematic uncertainties in A^max are dominated by the hadronic interaction models rather 
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FIG. 1: Distribution of A'max for 5,000 showers initiated by proton, iron nuclei and gammas at 
-|^gi9.5 Simulation has been done with the CONEX program. The thick solid line shows the 
sum of the A'max distributions for proton, iron nuclei and gamma primaries. 

than mass composition if proton and iron nuclei are considered, however the fluctuations 
due to both effects are very similar, 4.6% and 3.6%, respectively at 10^^'^ eV. If primary 
gamma are considered, mass composition would be more relevant than hadronic interaction 
models to the A^max fluctuations. 

In order to exemplify the increase in the uncertainty by allowing the introduction of 
photons in the cosmic ray flux, we also consider the peculiar case of equal abundances of 
protons, iron and gammas. Fig. IHl shows 5,000 air showers induced by gamma, proton and 
iron nuclei simulated with QGSJET and SIBYLL hadronic interaction models, at the energy 
of 10^^'^ eV. The total A'max uncertainty is ~ 7% which is comparable to the systematic 
uncertainties introduced by the missing energy correction calculation. 

III. SHOWER AND TELESCOPE SIMULATION 

The results presented in the current section have been obtained by using the CORSIKA 
6.2 simulation program. In addition, gamma showers have been simulated by considering 
the pre-shower effect. The thinning factor jl^ of the simulations is 10~® and the longitudinal 
air shower profiles were sampled in steps of 5 g/cm^. The energy thresholds were set to 0.1 
MeV, for electrons and photons, 0.3 GeV for hadrons and 0.7 GeV for muons. Protons and 
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FIG. 2: Distribution of A'^max for 5,000 showers simulated by QGSJET and SIBYLL. Simulation 
has been done with the CONEX program. The left hand side shows the primary protons at 10^^'^ 
eV while the right hand side illustrates the primary iron nuclei at the energy of 10^'^'^ eV. 

gamma ray showers were simulated in the energy interval from 10^^ to 10^°'^ eV in steps of 
0.1 dex. For each energy and primary particle we simulated 100 events which were recycled 
50 times by randomly drawing a new arrival direction and core position. For gamma showers 
simulated with the pre-shower effect only the core position was randomized since the relative 
direction of the primary particle with the geomagnetic field affects the development of the 
shower. 

The simulations of the HiRes telescopes were performed following the same general pro- 
cedure adopted by the HiRes Collaboration as published in Q| and Comparisons of 
our simulation with the HiRes Collaboration simulation and the HiRes data can be seen in 

A similar simulation program was written for the EUSO telescope in which we incorpo- 
rated the geometry details of that telescope. The EUSO experiment has been proposed to 
operate at 430 km of altitude looking downwards to measure the fluorescence light produced 
by the passage of a cosmic ray shower. It is configured with a 60 degree field of view covered 
by pixels of 0.1 degree resulting in 200.000 km^ detection area. 

For the HiRes studies we have simulated shower cores within a radius of 50 km from 
the telescope. The zenith angle has been randomly chosen from to 60 degrees. For the 
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FIG. 3: Distribution of A^max for 5,000 showers simulated by QGSJET and SIBYLL and initiated 
by proton, gamma and iron nuclei at 10^^'^ eV . Simulation has been done with the CONEX 
program. 

EUSO studies we allowed showers with cores within a radius of 250 km from the axis of the 
telescope and the zenith angle was chosen from to 90 degrees. 

Besides the configuration details, included in our simulation program, we also paid careful 
attention to the simulation of the atmosphere since the field of view of the EUSO telescope 
spans the entire atmosphere. In order to verify our code we have the EUSO acceptance as 
a function of the arrival direction of the primary particle. Fig. |3] illustrates the acceptance 
for the EUSO telescope, which is in good agreement with the predicted values given in P|. 

The importance of a detailed telescope simulation in the studies of A^max as an energy 
estimator can be seen in Fig. El in which the simulated A^max distribution as predicted by the 
CORSIKA code is compared to the reconstructed A^max, obtained after the detection by the 
EUSO and HiRes-I telescopes. It is evident from the figure that the distortion of the A^max 
distribution after the detection and reconstruction of the shower must be taken into account 
properly. Fluorescence telescopes have detection biases regarding arrival direction, primary 
energy and core position which are very difficult to control. Moreover, the reconstruction 
procedure also introduces fluctuations which depend on several parameters. These effects 
result in a distortion of the A^max distribution that depends on the configuration of the 
telescope and the chosen reconstruction method. Therefore the simulation of the telescope 
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FIG. 4: EUSO acceptance for proton showers at 10^'^ eV as a function of the zenith angle according 
to our simulation program. The heavy smooth line shows the isotropic distribution of events, the 
full line histogram corresponds to the events which trigger the telescope and the light line histogram 
the events which survived to the quality cuts. 
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FIG. 5: Comparison between the reconstructed and simulated A'max for the HiRes-I and EUSO 
telescopes. 

and the application of the reconstruction procedure is fundamental in order to obtain a 
meaningful relation between A^max and energy and to properly estimate the reconstruction 
errors attained by this method. 
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IV. PRIMARY ENERGY RECONSTRUCTION 



For the energy reconstruction of a shower there is one procedure which is unavoidable by 
all proposed methods, especially those in investigation here. Each pixel hit in the telescopes 
corresponds to a given atmospheric depth and the respective signal measured by it can be 
correlated to the number of particles or, more directly, to the energy deposited by the shower 
at that depth. From a limited number of points along the longitudinal profile, a function 
must be fitted in order to calculate, by extrapolation, the energy deposited by the shower 
beyond the limits of the field of view of the telescope. The most widely used function is the 
Gaisser-Hillas, which efficiency have been extensively tested and confirmed 

In this section, we are going to investigate the infiuence of the fitting procedure in the 
uncertainty of the energy for two reconstruction methods. The analysis will be done for a 
ground based experiment and may differ for a space based telescope. 

The standard method used by the HiRes and Auger Collaborations is based on the 
integration of the fitted Gaisser-Hillas function to obtain the shower energy and, from now 
on, it will be referred as the integral method. In addition, the method based on the A^max 
parameter will be referred as the A^max method. Our intention here is to evaluate the 
infiuence of the knowledge of a limited portion of the shower longitudinal development in the 
determination of the parameters related to energy, which are the integral of the longitudinal 
profile for the integral method and the A^max parameter for the iVmax method. 

We try to determine which parameter is less affected by the fact that the telescopes, in 
most of the cases, measure only a small part of the shower, and whether the extrapolation 
of the observed longitudinal profile depends on geometric parameters. In order to do so, 
we simulated 5,000 CORSIKA proton showers with energy 10^^'^ eV through the HiRes-II 
telescopes. For each shower we compare the integral of the longitudinal number of particles 
as calculated by CORSIKA using the whole profile of the shower with the integral of the 
fitted longitudinal profile based on the points actually measured by the telescope. The same 
comparison is made between the A^max given by CORSIKA based on the whole development 
of the shower and the reconstructed A^max after the restriction of the shower in the field 
of view of the telescope. Fig. IHl and [7| illustrate this procedure. These figures show the 
longitudinal profile simulated by CORSIKA in solid black lines. Star symbol shows the 
points measured by the telescope after the inclusion of noise while the dashed line shows the 
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Gaisser-Hillas fit to the points. 

The Gaisser-Hillas fit was implemented according to the procedures adopted by the HiRes 
Collaboration in 

and 0- 

Fig. ISl shows the distribution of errors in the determination of the integral of the lon- 
gitudinal profile and in the calculation of the A^max parameter for proton showers at 10^^'^ 
and lO^o-s eV. Table HI shows the mean and RMS values for all the distributions. The error 
is defined as the reconstructed value minus the simulated value divided by the simulated 
value. 

The same behavior is noticed in the A^max and integral errors. At 10^^'^ eV, the A^max 
parameter and the integral of the longitudinal profile are underestimated and the opposite 
happens at 10^°'^ eV where both values are overestimated. This might happen because show- 
ers with higher energy develop deeper in the atmosphere hitting ground level in an earlier 
stage. However, its is noticeable that the A^max parameter is less affected by reconstruction 
showing a mean closer to zero and a slightly smaller RMS value for both energies. 

The RMS values shown in Table IJ are the intrinsic fluctuations of the determination of 
the A^max parameter and of the integral value of the longitudinal profile. Nevertheless, we 
have not applied any quality cuts to the showers analyzed so far, besides asking that X^ax 
falls inside the field of view of the telescope. 

It has been claimed that A^max, as an energy estimator, should be superior to the integral 
method for showers with short track length detected by fluorescence telescopes. Other 
authors have also mentioned the dependence of the method on the shower zenith angle and 
on the distance from the telescope to the shower core. However, both zenith angle and 
distance to the core, are the variables which determine the track length of the shower in the 
field of view of the telescope. 

As illustrated in Figs. IHl and [7| the path length is the parameter which influences directly 
the fit of the Gaisser-Hillas function and, therefore, the energy reconstruction. The distance 
between telescope and shower also influences the energy reconstruction via the lateral size 
of the shower in the photomultiplier camera, as discussed in 0] and [2^, but this effect 
should be small compared to the reduction of the path seen by the telescope. 

Fig. El illustrates the reconstruction error of the integral and of the A^max parameter as a 
function of path length for proton showers with energy 10^^'^ eV. The line is a polynomial 
fit to the points. In this case, a slight better resolution for the A^max parameter is seen for 
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Integral 




Mean 


RMS 


Mean 


RMS 


1019.5 


0.02 


0.16 


0.03 


0.17 


1020.5 


-0.04 


0.10 


-0.06 


0.14 



TABLE I: Mean and RMS values for the error distributions given in Fig. |S1 




FIG. 6: Longitudinal development of a proton shower with energy 10^^'^ eV as given by CORSIKA 
(full line), as detected by the HiRes-II telescope according to our simulation where noise has been 
added (points) and as reconstructed by fitting a Gaisser-Hillas function (dotted line). 

the entire range of path length. 

Fig. shows the same plot for proton showers at 10^°'^ eV where a small difference 
between the reconstruction of the A^max and the integral can be seen for the range of path 
length from 5 to 35 degrees. 

These results show that the reconstruction of the Aniax parameter is equivalent to the 
reconstruction of the integral independently of the size of the shower detected by the tele- 
scope. Our calculations suggest that the reconstruction of the Amax parameter might become 
slightly more accurate than the integral reconstruction for the highest energies above ~ 10^° 
eV. 
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FIG. 7: Longitudinal development of a proton shower with energy 10^^'^ eV as given by CORSIKA 
(full line), as detected by the HiRes-II telescope according to our simulation where noise has been 
added (points) and as reconstructed by fitting a Gaisser-Hillas function (dotted line). 




FIG. 8: Error in the integral and A^max reconstructions due to the fit of a Gaisser-Hillas function 
using a limited number of points. 

V. RELATION BETWEEN PRIMARY ENERGY AND A^^ax 

Finally, we used the simulation programs already described to determine the relationship 
between the reconstructed A^max and the primary energy of the shower. We used the same 
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FIG. 9: Error in the integral and A'max reconstruction due to the fit of a Gaisser-Hillas function 
using a hmited number of points as a function of the path length for proton shower with energy 
-|^q19.5 Lines shows polynomial fit to the points in the range from 5 to 32 degrees. 

set of CORSIKA showers described above and simulated them through the HiRes-II and 
EUSO telescopes. The present are intended to assess the potential of the A^max parameter 
as an energy estimator by means of a full simulation and reconstruction chain including all 
types of systematic and statistical errors. 

The HiRes-II and EUSO telescopes were chosen to illustrate the method for ground and 
space based telescopes. The HiRes-II telescope is, in principle, very similar to the Auger 
telescopes regarding the variables explored here. Both experiments have 30 degrees of field 
of view in zenith and the pixel size for the Auger telescope is 1.5 degree while for the HiRes-II 
telescope is 1 degree. These small differences should not affect our conclusions. 



VI. HIRES TELESCOPE 



The reconstruction of the shower was also carried out according to the explanations 
given by the HiRes Collaboration in H and Q of their own procedures. The Xq and A 
parameters of the Gaisser-Hillas were fixed to -60 and 70 g/cm^ respectively. 

The following quality cuts were applied to the showers after reconstruction: 

• Angular speed < 11 



14 




FIG. 10: Error in the integral and A'max reconstruction due to the fit of a Gaisser-Hillas function 
using a limited number of points as a function of the path length for proton shower with energy 
-|^q19.5 Lines shows polynomial fit to the points in the range from 5 to 35 degrees. 

• Selected tubes < 7. 

• 0.85 < Tubes/degree < 3.0. 

• Photoelectrons/degree > 25. 

• Track length > 7°, or > 10° for events extending above 17° elevation. 

• Zenith angle < 60°. 

• 150 < Xinax < 1200 g/cm^, and is visible in detector 

• Geometry /n.d.f. < 10. 

• Profile fit x^/n.d.f. < 10. 

At each energy 5,000 showers were simulated and A^max reconstructed. Fig. ITTl illustrates 
the relation between the median A^max and the primary energy. The error bars represent the 
68% of confidence level of the reconstructed A^max distribution. The relationship found is 
almost linear and can be well fitted by the equation Logio(A'max) = — 8.75 + 0.97xLogio(-E). 

In order to test the reconstruction of the energy based on the A^max parameter we sim- 
ulated a second set of 5,000 independent showers. The energy was reconstructed using the 
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previous equation and by integrating the Gaisser-Hillas function. Fig. [represents the error 
in the energy reconstruction as a function of the primary energy. As suggested in the calcu- 
lations above, the A^max estimator is as good as the integral procedure for energies between 
10^^ and 10^^'^ eV. For energies above 10^^'^ eV the A^max estimator gets slightly superior to 
the standard integral method. 

We also studied the relationship between Amax and the primary energy for primary gamma 
showers. The relation is also linear and very similar to the one for protons shown in Fig. ITTl 
Fig. IT^ shows the relative difference between the median A^max for proton and gamma show- 
ers. The difference is defined as (A'max ~ ^max)/^max- The points in Fig. Uni show the 
difference calculated with the median A'max given by the reconstructed A'max distribution 
at each energy for proton and gamma primaries. The solid line is the difference calculated 
using the fitted relations between energy and A'max- 

For the energy range 10^^ ° to 10^°'^ eV, the difference is smaller than 12%. This difference 
would influence the energy reconstruction in a mixed composition of proton and gamma 
primaries. The same effect would be seen in both reconstruction methods: A^max and integral 
according to the discussion in Section |nl 

If a mixed composition of 50% proton and 50% gammas is considered the energy resolution 
would degrade but still would show a behavior similar to the one show in Fig. ^| At 10^^ 
eV the resolution would degrade to 20% for both methods. At lO^'^'^ the resolution for the 
A^max niethod would be 16% and for the integral procedure it would be around 18%. 

VII. EUSO 

The same analysis were extended to the EUSO telescope. The details of this experiment 
were included in the simulations as explained in Section lllll However, not only the tele- 
scope is different but its peculiar field of view also requires a considerable change in the 
shower reconstruction program. This can be verified in Fig. 0] in which we show the EUSO 
acceptance. 

Since there are no quality cuts yet optimized for the EUSO telescope, we imposed very 
loose ones requiring total path length greater than 0.6° and greater than 200 g/cm^ and 
Xmax in the field of view of the telescope. We would like to stress that different quality cuts 
will modify the relation between A'max and energy presented here. 
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FIG. 11: Relation between the reconstructed A^max and the primary energy for showers initiated 
by protons. Results shown for the HiRes-II telescope. 
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FIG. 12: Error in the energy reconstruction by the integral and A^max methods as a function of the 
primary energy. Results shown for the HiRes-II telescope. 

Fig. HH shows the relation between the reconstructed A'max and the simulated energy, 
which can be well approximated by the equation: Logio(A^max) = -8.95 + 0.90 x Logio(-E'). 

In order to test the reconstruction of the energy based on A^max we simulated a second set 
of 5,000 independent shower. The energy was reconstructed using both, the equation above 
and the integration of the Gaisser-Hillas fitted profile. Fig. ITHl shows the error in the energy 
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FIG. 13: Relative difference between the median N^nax for proton and gamma showers. Points 
shows the difference calculated with the median of the reconstructed Nj^^ distribution and the 
solid line shows the difference calculated with the fitted relation between energy and Nmax- Note 
that the line is not a linear fit to the points. 

reconstruction as a function of the primary energy. For energies below 10^°-^ eV, the A^max 
estimator was superior to the integral method. For energies above 10^°'^ eV, both methods 
are equivalent. 

These results are different from those for the HiRes-II telescopes, which confirms that 
the relation between iVmax and primary energy depends on the telescopes and details of the 
analysis, mainly the quality cuts. 

The exact influence of each quality cut on the energy resolution is not completely un- 
derstood at present and is out of the scope of this paper. However, the results shown here 
suggest that the quality cuts imposed by the HiRes Collaboration are well optimized in order 
to guarantee a good energy resolution for the integral method for the energy ranging from 
10^^ to 10^°'^ eV. On the other hand, the loose cuts applied by us to the EUSO analysis do 
not achieve an optimal selection of event for energies below 10^°'^ eV resulting in an increase 
of the reconstruction error for both techniques. It is also apparent that A^max, for EUSO, is 
less sensitive to well cahbrated cuts because, despite the fact that the errors in this method 
increases at lower energies, they still remain within reasonable bounds, i.e., smaller than 
15% down to 10^9 eV. 
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FIG. 14: Relation between the reconstructed A^max and the simulated primary energy for showers 
initiated by protons. Results are for the EUSO telescope. 




FIG. 15: Error in the energy reconstruction by the integral and A^max methods as a function of the 
primary energy. Results are for the EUSO telescope. 

VIII. CONCLUSION 

In the present work, we tested the number of particles at the maximum air shower de- 
velopment, A^max, as an energy estimator. The uncertainties associated with the simulation 
models were calculated and compared to the uncertainties in the missing energy calculation. 
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A similar intrinsic fluctuation of around 5-10% at 10^^'^ eV was verifled for both approaches. 
These uncertainties decrease with increasing energy. 

The systematic uncertainties in the determination of A^max are dominated by mass compo- 
sition rather than hadronic interaction models if gamma showers are taken into account. If 
gamma showers are excluded from the analysis, mass composition and hadronic interaction 
models would contribute equivalently to the total uncertainties. 

Full telescope simulations were considered in our analysis. The flrst source of uncertainty 
in the reconstruction procedure is the flt of the points measured by the telescope using a 
Gaisser-Hillas function. We investigated the statistical uncertainties involved in the eval- 
uation of A/'max and of the integral of the longitudinal distribution function as a function 
of the path length seen by the telescope. A'max appears to be slightly better reconstructed 
independent from the path length observed by the telescope. At 10^^'^ eV the difference in 
reconstructing A^max and the integral is very small, however with increasing energy, at 10^°'^ 
eV for example, the difference is of the order of a few percent, depending on energy (see 

Fig. cm). 

The relation between the A^^ax and primary energy for proton showers was calculated 
for the HiRes-II telescope including the quality cuts used by the HiRes collaboration. A^max 
appears to be equivalent to the standard integral procedure as an energy estimator, with 
a small superiority of the A^max approach for the highest energies (see Fig. IT^ . We also 
calculated the same relation for gamma showers and the relative difference for the proton 
shower is smaller than 12%. If a mixed composition of 50% proton and 50% gammas is 
considered the energy resolution degrades to 20% at 10^^ eV for both methods. At 10^°'^ 
the resolution for the A'max method is 16% while for the integral procedure is ~ 18%. 

The same analysis were performed for the EUSO telescope and we demonstrated that the 
A'max method also works well for space based telescopes. The speciflc relationships, however, 
depend on the telescopes characteristics and reconstruction methods. 

Finally, A'max was established as a good energy estimator for fluorescence telescope with 
some advantages over the standard integral procedure. 
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